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Abstract 

In this paper we investigate the role of Parodi's relation in the well-posedness and sta- 
bility of the general Ericksen-Leslie system modeling nematic liquid crystal flows. First, 
we give a formal physical derivation of the Ericksen-Leslie system through an appropriate 
energy variational approach under Parodi's relation, in which we can distinguish the conser- 
vative/dissipative parts of the induced elastic stress. Next, we prove global well-posedness 
and long-time behavior of the Ericksen-Leslie system under the assumption that the viscosity 
P4 is sufficiently large. Finally, under Parodi's relation, we show the global well-posedness 
and Lyapunov stability for the Ericksen-Leslie system near local energy minimizers. The 
connection between Parodi's relation and linear stability of the Ericksen-Leslie system is 
also discussed. 

Keywords: Liquid crystal flows, Ericksen-Leslie System, Parodi's relation, uniqueness of 
asymptotic limit, stability. 
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1 Introduction 

Liquid crystal is often viewed as the fourth state of the matter besides the gas, liquid and solid, 
or as an intermediate state between liquid and solid. It possesses none or partial positional order 
but displays an orientational order at the same time. The nematic phase is the simplest among 
all liquid crystal phases and is close to the liquid phase. The molecules float around as in a liquid 
phase, but have the tendency of aligning along a preferred direction due to their orientation. 
The hydrodynamic theory of liquid crystals due to Ericken and Leslie was developed around 
1960's [12, 13, 25, 26j. Earlier attempts on rigorous mathematical analysis of the Ericksen-Leslie 
system were made recently |32j (see [28 , 30, 31J for a simplified system which carried important 
mathematical difficulties of the original Ericksen-Leslie system, except the kinematic transport 
of the director field). 

The full Ericksen-Leslie system consists of the following equations (cf. [14, 26j[2Tll32]): 

Pt + vVp = 0, (1.1) 
pv = P F + V • a, (1.2) 
piu = piG + g + V -7T. (1.3) 
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Equations (jl.lf ) f|1.3|) represent the conservation of mass, linear momentum, and angular mo- 
mentum, respectively, with the anisotropic feature of liquid crystal materials exhibited in (|1.3|) 
and its nonlinear coupling in (|1.2p (cf. [26,32]). In this paper, we consider the flow of an in- 
compressible material, namely, V • v = 0. Here, p is the fluid density, p\ is a (positive) inertial 
constant, v = (vi,V2,vs) T is the flow velocity, d = (d\,d2,ds) T is the orientational order pa- 
rameter representing the macroscopic average of the molecular directors, g is the intrinsic force 
associated with d, ir is the director stress, F and G are external body force and external director 
body force, respectively. The superposed dot denotes the material derivative dt + v ■ V. The 
notations 

A=^(Vv + V T v), n = ^(Vv-V T v), 

ui = d = d t + (v ■ V)d, N = u-Qd, 

represent the rate of strain tensor, skew-symmetric part of the strain rate, the material derivative 
of d (transport of center of mass) and rigid rotation part of director changing rate by fluid 
vorticity, respectively. 

We have the following constitutive relations in the system (]1.1|) - (]1.3|) for a, tt and g: 

dW 

&ij = -PSij - p-^i—dk,j + (Tij, (1.4) 
dW 

nj = fcdj + p—, (1.5) 

dW 

9i = idi - fijdij - p— + gi. (1.6) 

Here P is a scalar function representing the pressure. The vector /3 = (/3i, /?2, Pz) T an d the scalar 
function 7 (sometimes called director tension) are Lagrangian multipliers for the constraint on 
the length of director such that \d\ = 1, with the Oseen-Frank energy functional W for the 
equilibrium configuration of a unit director field: 

W = ^(V-d) 2 + ^\dx(Vxd)\ 2 + ^-\d-(Vxd)\ 2 

+(k 2 + k 4 )[tr(Vd) 2 - (V-d) 2 }. (1.7) 

We note that the forth term in (|1.7p 

tr(Vd) 2 - (V • df = V • [(Vd)d - (V • d)d] 

is a null Lagrangian, which only depends on the value of the trace of d on the boundary (cf. PQ). 
The kinematic transport of the director d (denoted by g) is given by: 

gi = XiNi + X 2 djAji = Ai (V, + ^djAj^J , (1.8) 

which represents the effect of macroscopic flow field on the microscopic structure. The material 
coefficients Ai and A2 reflects the molecular shape (Jeffrey's orbit [21]) and the slippery between 
the fluid and the particles (see discussions in Section [3|) . The first term of (jl.8p represents the 
rigid rotation of the molecule, while the second term stands for the stretching of the molecule 
by the flow. 

The stress tensor a has the following form: 



= pidkAkpdpdidj + p 2 Nidj + p^Nj + p 4 A. 
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(1.9) 



These (independent) coefficients /Ui, which may depend on material and temperature, are 

usually called Leslie coefficients. These coefficients are related to certain local correlations in 
the fluid (cf. |llj). For convenience, /i^s are called viscous coefficients in later sections. 

In order to handle the higher-order nonlinearities due to the nonlinear constraint \d\ = 1 (i.e., 
the Lagrangian multipliers /?, 7), one can introduce a penalty (or relaxation) approximation of 
Ginzbug-Landau type, by adding one term 

Hd) = ~^{\d\ 2 ~ 1? 

in W. Physically this term can be attributed to the extensibility of the molecules. After the 
discussions for each e > 0, we then take the limit as e — > 0. This method is motivated by the work 
on the gradient flow of harmonic maps into the sphere (see, e.g., [S]), but whether the solution 
of the Ericksen-Leslie system with Ginzburg-Landau approximation converges to that of the 
original one with constraint \d\ = 1 as e tends to zero is still a challenging problem. Nevertheless, 
the reformulated system with penalty approximation also has natural physical interpretations. 
It is similar to what Leslie proposed in |27j for the flow of an anisotropic liquid with varying 
director length. Mathematically, it can also be related to models for nematic liquid crystals 
with variable degree of orientation proposed by Ericksen in [14J under specific conditions. In 
particular, {x : d(x, t) = 0} represents the transition region of isotropic fluids. We refer to [30] 
for more discussions. 

For simplicity, in this paper, we focus on the relaxation form of the elastic energy associated 
with d: 

W{d)= l -\Vd\ 2 + ^{\d\ 2 -lf. (1.10) 

It is obvious that this choice of W corresponds to the elastically isotropic situation, i.e., taking 
k\ = &2 = £3 = 1, &4 = in (|1.7j) . The corresponding problem with general Oseen-Frank energy 
(|1.7p can be treated in a similar way, but the argument is more involved. With the choice of the 
penalized energy W, we can remove the Lagrangian multipliers and set 7 = f3j = 0. Since the 
inertial constant p\ is usually very small (cf. [10J ) , we take p\ = 0. Moreover, we assume that 
the density is constant (which in turn yields the incompressibility V • v = 0) and the external 
forces vanish, i.e., p = 1, F = 0, G = (cf. [32]). Note that in the incompressible cases the 
assumption F = means that there is no exterior nonconservative forces. 
Thus, the full Ericksen-Leslie system (|l.ip - (|1.3p can be reformulated to 

v t + vVv + VP = -V • (VdOVd) + V • a, (1.11) 
V-u = 0, (1.12) 

d t + (v-V)d-nd+^Ad = (Ad- /(d)), (1.13) 

Ai Ai 



where 
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f{d) = F'{d) = - 2 {\d\ 2 -l)d 

and a is given by (|1.9|) . We denote by VdQVd the 3 x 3-matrix whose (i, j)-entry is Vid ■ Vjd, 
1 < i,j < 3. In the following text, we just set e = 1 and our results indeed hold for any arbitrary 
but fixed e > 0. In this paper, we will focus on the bulk properties of the Ericksen-Leslie system. 
For this, we consider the equations (|l.lip - (|1.13p subject to periodic boundary conditions (i.e., 
in torus T 3 ): 

v(x + ei,t) = v(x,t), d(x + ei,t) = d(x,t), for (x, t) € dQ x M + , (1-14) 
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and initial conditions 



v\t=o = vq(x), with V • v = 0, d\ t =o = do(x), for x G Q 



(1.15) 



where Q is a unit square in M 3 . 

Due to temperature dependence of the Leslie coefficients, there exists different behavior 
between various coefficients (cf. [H]): Ha- which does not involve the alignment properties- is a 
rather smooth function of temperature; but all the other fjfs describe couplings between the 
molecule orientation and the flow, and are thus affected by a decrease in the nematic order \d\. 
In this paper, we just look at the isothermal case where fi's are assumed to be constants. The 
following relations are frequently introduced in the literature (cf. [26^127]) 



Relations given in (|1.16|) are necessary conditions in order to satisfy the equation of motion 
identically (cf. [261 Section 6]). f)1.17|) is called Parodi's relation (cf. [10]), which is derived 
from Onsager reciprocal relations expressing the equality of certain relations between flows and 
forces in thermodynamic systems out of equilibrium (cf. |37|). Under the assumption of Parodi's 
relation, we see that the dynamics of an incompressible nematic liquid crystal flow involve five 
independent Leslie coefficients in (|1.9p . 

Since the mathematical structure of the Ericksen-Leslie system is quite complicated, past 
existing work was almost all restricted to its simplified versions (cf. [31 [301 [103 [MIES EI]). As 
far as the general Ericksen-Leslie system is concerned, there is few known result in analysis (cf. 
e.g., [32J). In [32J, well-posedness of the general Ericksen-Leslie system (|l.lip - ([1.13p subject to 
Dirichlet boundary conditions was proved under the special assumption A2 = 0, which imposed 
an extra constraint on those Leslie coefficients. This physically indicates that the stretching 
due to the flow field is neglected, which is more feasible for small molecules. Mathematically 
this assumption brings great convenience since a weak maximum principle for \d\ holds (cf. [321 
Theorem 3.1]). For the general system (jl.lip - (|1.13p . the maximum principle for \d\ fails when 
A2 7^ 0. This leads to extra difficulties in the study of well-posedness, especially in dealing with 
those highly nonlinear stress terms in a (cf. [321144] ). Even in the 2D case, it is hard to obtain 
global existence of solutions without any further restriction on these viscous coefficients. This 
is rather different from regular Newtonian fluid cases. 

Summary of results. The purpose of this paper is to study the connections between physical 
parameters, namely, the Leslie coefficients and the well-posedness as well as stability proper- 
ties of the general Ericksen-Leslie system (jl.lip - (|1.13p . In particular, we focus on the role of 
Parodi's relation ()1.17p . Parodi's relation is a consequence of Onsager 's reciprocal relations in 
the microscale descriptions of liquid crystals [371[38] , which are nevertheless independent of the 
second law of thermodynamics. Although the physical interpretation of the reciprocal relation 
is related to microscopic reversibility and laws of detailed balance of microscopic dynamics [36], 
the thermodynamic basis of Onsager's reciprocal relations have been discussed and debated by 
many researchers (cf. [46 J. There arc evidences that for particular materials, Onsager's relations 
and their counterparts may serve as stability conditions (cf. OdS]). In this paper, we provide 
specific mathematical verifications for the nematic liquid crystal flow. 

First, in Theorems 14.11 and 14.21 under the assumption that the fluid viscosity ^4 is sufficiently 
large, we show existence and uniqueness of global solutions within suitable regularity classes 
and their long-time behavior (uniqueness of asymptotic limit). In this case, we see that Parodi's 
relation is not a necessary assumption for the well-posedness and long-time dynamics, while the 
large viscosity constant [it± plays a dominative role. 

Next, without the largeness assumption on we first prove local wellposedness of the 
Ericksen-Leslie system (|l.lip - (|1.13p (cf. Theorem 15. ip . Furthermore, in Theorem 15.21 we prove 



Ai = /U 2 - fJ-3, A 2 = - /i 6 
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global well-posedness and Lyapunov stability of the Ericksen-Leslie system, when the initial data 
is near certain equilibrium (local minimizer of the elastic energy W given by (jl.iup ). We see that 
Parodi's relation turns out to be crucial (as a sufficient condition) in obtaining well-posedness 
and (nonlinear) stability of the Ericksen-Leslie system. 

Finally, we demonstrate the connection between Parodi's relation and linear stability of the 
original Ericksen-Leslie system (|l.ip - (|1.3p (namely, with the constraint |d| = 1). The result 
obtained in Theorem 16.11 indicates that without Parodi's relation, the linearized Ericksen-Leslie 
system admits unstable plane wave solutions. In other words, Parodi's relation is a necessary 
condition for linear stability of the Ericksen-Leslie system. 

Remark 1.1. Our results presented in this paper are stated in the three dimensional case n = 3. 
When the spatial dimension n = 2, if we consider the velocity field v : Q x [0, T] — > M? and 
director d : Q x [0, T] — > M 2 , namely, the molecule director d is also confined in a plane, then it 
is easy to verify that all the results we obtained below for the 3D case also hold in 2D (sometimes 
even under weaker assumptions, see e.g., Remark |^.i| ). However, if one wishes to consider the 
Ericksen-Leslie system in a 2D domian Q, C M 2 but the director field d is still allowed to be a 
three dimensional vector, some troubles will come up. For instance, the parallel transport terms 
£ld (rotation) and Ad (stretching) cannot be properly defined, because f2 and A are 2x2 matrices 
but d is a 3D vector. We want to mention that such problem does not apply to simplified liquid 
crystal system of small molecules }30f . In particular, we refer to recent works f29, 33, 43jj for 
a simplified liquid crystal system in 2D but the director d : x [0, T] — > S 2 , which is three 
dimensional (with the constraint \d\ = 1). 

Plan of the paper. The remaining part of the paper is organized as follows. In Section 2, 
we discuss the basic energy dissipation law of the system (|l.lip - (|1.13p . In Section 3, for the 
given energy law, we re-derive the Ericksen-Leslie system via an energy variational approach. 
In particular, under Parodi's relation, we show the specific relations between results from Least 
Action Principle and those from Maximum Dissipation Principle. In Section 4, we prove global 
well-posedness under large viscosity assumption on n± and the long-time behavior of global 
solutions. In particular, we show that any global solution will converge to a single steady state 
as time tends to infinity and provide an estimate on the convergence rate. In Section 5, we prove 
the well-posedness and stability when the initial velocity is near zero and the initial director is 
close to a local energy minimizer under Parodi's relation. In Section 6, we discuss the connection 
between Parodi's relation and linear stability of the original Ericksen-Leslie system. In Section 
7, the appendix section, we present some detailed calculations needed for the previous sections. 

2 Basic Energy Law 

Generally speaking, singularities that can be observed for a physical system are those ener- 
getically admissible ones (cf. [31]) • It has been pointed out that the Ericksen-Leslie system 
(|l.lip - ()1.15p obeys some dissipative energy inequality under proper assumptions on those phys- 
ical coefficients (cf. [52"]). 

The total energy of the Ericksen-Leslie system (jl.lip - (|1.15p consists of kinetic and potential 
energies and it is given by 



For the sake of simplicity, we denote the inner product on L 2 (Q) (or L 2 (Q), for the corresponding 
vector space) by (•, •) and the associated norm by || • ||. 




Q 



(2.1) 
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By a direct calculation with smooth solutions (v, d) to the system ([l.lip - ([1.15p . we have 
(cf. [321 Theorem 2.1] for detailed calculations for the corresponding initial boundary value 
problem) 

j£{t) = - [ [fn\d T Ad\ 2 + ^\Vv\ 2 + (^ + ^)\Ad\ 2 )dx 

+\!\\N\\ 2 + (A 2 - M2 - VtoKN, Ad). (2.2) 
Here and after, we always assume that 

Ax < 0, (2.3) 

+ > 0, (2.4) 

Mi > 0, fi 4 > 0. (2.5) 

These assumptions are assumed to provide necessary conditions for the dissipation of the di- 
rector field [15l[27j. As indicated in [32], the assumption (|1.16p guarantees the existence of the 
Lyapunov-type functional. However, we note that Parodi's relation (|1.17p is not necessary in 
the derivation of (|2.2p . If (|1.17j) is employed, i.e., A2 = —(fJ-2 + 1^3), we immediately arrive at 
the energy inequality obtained in [32 \ Theorem 2.1]. Moreover, if we further assume A2 = 0, it 
follows from (|2.2p - (|2.5p that £{t) is decreasing in time, which is exactly the case studied in 



Lemma 2.1 (Basic energy law with Parodi's relation). Suppose that the assumptions (|1.16p . 
(|1.17|) . (|2.3|) . (|2.4|) and (|2.5p are satisfied. In addition, if we assume 

^<M5+M6, (2.6) 
— M 

then the total energy £(t) is decreasing in time such that 

j t £{t) = -J^^Adl' + ^lVv^dx + j-WAd-fid)]] 2 

^+^+ { -^-)\\Ad\\ 2 
< 0. (2.7) 

Proof. By Parodi's relation (|1.17p . i.e., A2 = —{^2 + ^3)? we infer from the transport equation 
of d (cf. ([TT3]n that 



Ai||iV|| 2 + (A 2 -/i 2 -M3)(A r ,^) 
(N, AiiV + \ 2 Ad) + X 2 (N, Ad) 

(N, \\N + X 2 Ad) +(n + ^-Ad, X 2 Ad\ - ^-\\Ad\\' 



Ai 



N+^Ad 
Ai 



{X2)2 -\\Ad\\ 2 



A 



°P h\Ad- m f-^f\\Adf. (2.8) 

Inserting the above result into (|2.2p . we arrive at our conclusion. □ □ 

On the contrary, if Parodi's relation (|1.17p does not hold, additional assumptions have to be 
imposed to ensure the dissipation of the total energy. 
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Lemma 2.2 (Basic energy law without Parodi's relation). Suppose that (|1.16[) . (|2.3p . (j2.4[) and 
(|2.5|) are satisfied. If we further assume that 



|A 2 - - ^3 1 < 2 v / -AiVAi5 + M6, (2-9) 
iaen iae following energy inequality holds: 

-£{t) < - [ (m\d T Ad\ 2 + ^\Vv\ 2 )dx<0. (2.10) 



Moreover, if 



|A 2 - - ^3 1 < 2y / -Ai^//i 5 + /i 6 , (2.11) 



i/ien i/ie dissipation in (|2.10|) 6e stronger in the sense that there exists a small constant 
V > o, 



< 0. (2.12) 
Proof. The conclusion easily follows from (j2.2|) and the Cauchy-Schwarz inequality. □ □ 



3 Energy Variational Approaches 

The energy variational approaches (EnVarA) provide unified variational frameworks in studying 
complex fluids with microstructures (cf. [20]). From the energetic point of view, the Ericksen- 
Leslie system (|l,lip - ()1.13p exhibits competition between the macroscopic flow field and the 
microscopic director field, through the coupling between the kinematic transport of the director 
d by the macroscopic velocity field v and the averaged microscopic effect in the form of induced 
macroscopic elastic stress on the macroscopic flow field. This contributes to some interesting 
hydrodynamic and rheological properties of the liquid crystal flows. Based on the basic energy 
law in Section 2, and due to the special feature of nematic liquid crystal flow such that the 
molecular orientations are transported and deformed by the flow under parallel transport, we 
shall develop a formal physical derivation of the induced elastic stress through EnVarA. This 
will provide us with a further understanding of the competition between hydrodynamic kinetic 
energy and internal elastic energy due to the presence of the orientation field d. 

The energetic variational treatment of complex fluids starts with the energy dissipative law 
for the whole coupled system [2TJ1I5U]: 

dE tot 

— = - v ' 

where E tot = E kmettc + E tnt is the total energy consisting of the kinetic energy and free energy. 
Here T> is the dissipation function which is equal to the entropy production of the system in 
isothermal situations. Following Onsager's linear response assumption, we assume that T> is a 
linear combination of the squares of various rate functions such as velocity, rate of strain or the 
material derivative of internal variables (cf. [20|l37l - f39|l50| ) . The EnVarA combines the maximum 
dissipation principle (for long time dynamics) and the least action principle, or equivalently, the 
principle of virtual work (for intrinsic and short time dynamics) into a force balance law that 
expands the conservation law of momentum to include dissipation (cf. |6,23j). The least action 
principle gives us the Hamiltonian (reversible) part of the system related to conservative forces. 
Meanwhile, the maximum dissipation principle provides the dissipative (irreversible) part of 
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the system related to dissipative forces. In this way, we can distinguish the conservative and 
dissipative parts among the induced stress terms. 

In the context of basic mechanics, both hydrodynamics and elasticity, the basic variable is 
the flow map x(X, t) (particle trajectory for any fixed X) . Here, X is the original labeling 
(the Lagrangian coordinate) of the particle, which is also referred to as the material coordinate, 
while x is the current (Eulerian) coordinate and is also called the reference coordinate. For a 
given velocity field v(x,t), the flow map is defined by the ordinary differential equations: 

x t = v(x(X,t),t), x(X,0) = X. 

The deformation tensor F associated with the flow field is given by 

dxi 



dX 3 



Without ambiguity, we define ¥(x(X,t),t) = ¥(X,t). Applying the chain rule, we can see that 
¥(x,t) and ¥~ T (x,t) satisfy the following transport equations (cf. e.g., [T71l2"4"]) 

¥ t + v ■ VF = X7v¥, 

¥t T + v ■ V¥~ T = -V T v¥- T . 

Kinematic transport of the director field d represents the (microscopic) molecules moving in the 
(macroscopic) flow [24,44J. For general ellipsoid shaped liquid crystal molecules, the transport 
of d can be represented by 

d(x{X,t),t) =Ed (X) (3.1) 

with do(X) being the initial configuration. The deformation tensor ¥(x(X,t),t) carries all the 
information of micro structures and configurations. It satisfies the following transport equation 
whose right-hand side can also be reformulated into a combination of a symmetric part and a 
skew part: (cf. [SUSHIS]) 



E t + v ■ VE 



E 



aVv + (1 - a)(-V J v) 
fiE + (2a - 1)AE. (3.2) 



Such solutions are called Jeffrey's orbits (cf. |21j). By the fundamental work of Jeffrey [21], the 
parameter 

I? = 2a-l = ^G[-l,l], r€R 
r z + 1 

is related to the aspect ratio of the ellipsoids. Recently, we have shown that tj can also be related 
to the slippage between the particles and the flow |44j . In the present case, we see that 

1 / A 2 

In what follows, we shall apply EnVarA to recover the system (|l.lip - ([1.13p from the basic 
energy law under the assumption that both (|1.16p and Parodi's relation (|1.17|) are satisfied. The 
kinetic energy E kinetlc and internal elastic energy E int of the system (jl.lip - (|1.13p are given by 

E kinetic = -\\v\\ 2 , E int = E(d) = -\\Vd\\ 2 + [ T(d)dx. 
2 2 J Q 

The Legendre transformation yields the action functional A of the particle trajectories in terms 
of the flow map x(X, t): 



A(x) = f (E kmetic - E int )dt, 
Jo 
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which represents the competition between the kinetic energy and the internal energy. If there is 
no internal microscopic damping, we deduce the total (pure) transport equation of d from 
such that 



Dd 

151 



dt + v-Vd- aVv d + (1 - a)(V T v)d 
A 2 

ck + v ■ Vd - 0.d + -^Ad 
Ai 



(3.3) 



The least action principle optimizes the action A with respect to all volume preserving trajec- 
tories x(X,t), i.e., 5 X A = 0, with incompressibility of the fluid V • v = 0. Consequently, we 
obtain the conservative force balance equation of classical Hamiltonian mechanics (see (|7.6p for 
its weak variational form) 



Vt + vVv = -VP - V • (Vd Vd) + V • or, 



(3.4) 



where 



a 



I. ( x A 2 



A J (Ad-/(d))0d 

\2 



+ 2 [l + fJd^iAd-fid)). 



(3.5) 



Here, the symbol denotes the usual Kronecker multiplication, namely, (a <g> b)ij = aibj for 
a, b € M 3 and 1 < i,j < 3. We also note that the stress tensor cr is not symmetric due to the 
different coefficients of its two components. Together with (|3.3p . we recover the conservative 
(Hamiltonian) part of the full system fjl . 1 1 j) — (| 1 . 1 3j) (see Section [7711 for the detailed calculations). 

On the other hand, taking the internal dissipation into account together with the transport 
equation (j3.3|) . we get 



A 2 

dt + v- Vd-ttd+ -±Ad 
Ai 



hit 



1 5E 
Ai 6d 

-i(Ad 

Al 



f(d)), 



(3.6) 



which reflects the elastic relaxation dynamics. The dissipation functional T> to the system (jl.lip - 
(|1.13p is in terms of the variables A and N (cf. (|2.2p ) (we remark that our dissipation functional, 
like in [2j, departs from those loosely defined by Onsager in [3S])- Under Parodi's relation (|1.17|) 
and by (|3.6p . it can be transformed into the following form (cf. (|2.8p ) 



V 



^Adf + ^-WVvf- 
+ ^B + ^ + ^-)\\Ad\\ 



A 2 

d t + v ■ Vd - n d + -^A d 
Ai 



(3.7) 



According to the maximum dissipation principle [37H39] . we take S v (|^) = (performing 
variation with respect to the rate function, i.e., the velocity v in Eulerian coordinate) with 
incompressibility of the fluid V ■ v = 0. This yields the dissipative force balance law equivalent 
to the conservation of momentum (see Section 17.21 for the detailed calculations) : 







-VP - V • (Vd Vd) + V • a, 



(3.8) 
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where 

a = ni(d T Ad)d <g> d + fj, 2 N ® d + fi 3 d ® iV + [i^A 

+[i 5 Ad® d + [x 6 d® Ad, (3.9) 

with constants ^ 

M2=2( A i- A 2)> = -^ Xl + A2 )- 

Combining (|3.8p with (|3.6|) . we recover the dissipative part of the full system (|l.ll|) - (|1.13p . 
which stands for the macroscopic long time dynamics. 

The Ericksen-Leslie system (|l.lip - (|1.13p is the hybrid of these two conservative/dissipative 
systems. Combining the dissipative part derived from maximum dissipation principle (cf. (|3.8|) ) 
with the conservative part derived from the least action principle (cf. (|3.4p ). and taking into 
account the total equation of d (|3.6p . we recover the full system (jl.lip - (|1.13p . 

Remark 3.1. We first observe from ([33]) and (i.e., -X\N - X 2 Ad = Ad - f{d)) that 

a = [i 2 N ®d + n-sd <g> N + rj 5 A d®d + rj 6 d® Ad, (3.10) 



with constants 



M2 = ^(Ai-A 2 ), Hz = --(A x + A 2 ) 

1 (s _ 1 A , (A 



^2 



" 5 =2( A2 ^J- " e= -H A2+ "vJ- <3 - n) 

The interesting fact from the above derivation is that the induced stress terms 

-Vd @Vd + /j, 2 N ®d+ fi 3 d ® N + r] 5 Ad ®d + rj 6 d® Ad 

can be derived either by the least action principle ( cf. ()3.10|) ) or the maximum dissipation prin- 
ciple (contained in (I3.9p ). Therefore, they can either be recognized as conservative or dissipative. 
However, the remaining part in (|3,9p 



Hi{d T Ad)d ®d + mA + (fi 5 - n 5 )Ad <g> d + {(x & - rj 6 )d <g> Ad (3.12) 

can only be derived by the maximum dissipation principle. This fact indicates that these terms 
in (|3.12p are dissipative. In particular, at the critical value of \ 2 , i.e., 



|A 2 | = V-Aiv 7 ^ + M6, (3.13) 
the dissipation functional of the system ([l.lip (j!.13[) is reduced to 



V = ^\\d T Adf + ^\\Vv\\ 2 -X 1 



X 2 
d t + v - Vd-£ld+ -^A d 
Ai 



It turns out that /X5 = 775, [xq = tjq and the only two dissipative terms are given by those associated 
with fj,% and □ 

Finally, we look at some special cases of the system (|l.lip - (|1.13p . We assume that (|1,16D 
(|1.17p are satisfied and set 

Mi = 0, 

IJ-2 = ^ (Xi - X 2 ) , V3 = --{X! + X 2 ) , 
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1 /\ (A 2 ) 2 \ 1 /. (A 2 ) 2 

Since ()3.13j) is now satisfied, then the system f)l . 1 1 p — (|1 . 13[) can be reduced to 

v t + v ■ Vv + Vp = y At> - V • (Vd Vd) + V • a, (3.14) 

V-« = 0, (3.15) 

dt + v-Vd- ^-Vv d - ^-V T v d = - — (Ad - /), (3.16) 
Ai Ai Ai 

where 

a = -^-{Ad- f)®d - ^-d®(Ad- f). (3.17) 
Ai Ai 

Remark 3.2. The system (|3.14|) - (|3.17p is consistent with these simplified models studied in 

(1) The rod-like molecule model: 

M2 = Ai = -A 2 , = 0. 
In this case, the director field d satisfies the kinematic transport relation 

d(x(X,t),t) = Fd (X), where F = VvF. 

(2) The disc-like molecule model: 

(J>2 = 0, fi 3 = -Ai = -A 2 . 

In f/tis case, d satisfies 

d(x(X, t),t) = ¥~ T d (X), where ¥~ T = -V T v¥- T . 

(3) The sphere-like molecule model: 

M2 = y, M3 = -y, A 2 = 0. 

In i/iis case, d satisfies 

d(x(X,t),t) = Ed (X), where E = ^(Vu - V T v)E. 

4 Well-posedness and Long-time Behavior for Large Viscosity 

For any Banach space X, we denote by X the vector space (X) r , r £ N, endowed with the product 
norms. We recall the well established functional settings for periodic problems (cf. [45] ) : 





= {u 6 H m (M?, R) u(x + ei) = u(x)}, 




K(Q) 


= H™{Q) n iu : y u(x)dx = 1 , 




H 


= {f £ L 2 (Q), V • v = 0}, where L 2 (Q) 


= Hj(Q) 


V 


= {«€H*(Q), V •<; = ()}, 




V 


= the dual space of V. 
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We denote the inner product on Lp(Q) (or Lp(Q)) as well as H by (•, •) and the associated norm 
by || • || . The space H™(Q) will be short-handed by H™ and the H m -mner product (m G N) 
can be given by (v, u)u m = Y^\=o(D K v, D K u), where k = («i, n n ) is a multi-index of length 
\k\ = Y17=l K i an d D K = d£} , <9£™ . We denote by C the genetic constant possibly depending on 
A^s, n^s, Q, f and the initial data. Special dependence will be pointed out explicitly if necessary. 
Throughout the paper, the Einstein summation convention will be used. 

As mentioned in the introduction, we use the Ginzburg-Landau approximation to reduce 
the order of nonlinearities caused by the constraint \d\ = 1. We note that either for a highly 
simplified liquid crystal model (cf. [30]), or for the general Ericksen-Leslie system fll.ll)) — (|1.15|) 
with the artificial assumption A2 = (cf. |32j). a certain type of maximum principle holds for 
the d-equation, namely, if | do I < 1 then \d\ < 1. This fact still holds for our current periodic 
settings with the same assumption on A2. Then combing the basic energy law, one can deduce 
that 

v G L°°(0,T; H) D L 2 (0,T; V), (4.1) 
dG L°°(0,r;H^nL^)nL 2 (0,r;H2), (4.2) 

which is sufficient for the following formulation of weak solutions: 

Definition 4.1. (v, d) is called a weak solution of (jl.lip - (|1.15p in Qt = Q x (0, T) if it satisfies 
(|4.ip . (|4.2p and for any smooth function ifi(t) with ip(T) = and <p(x) 6 H^, the following weak 
formulation together with the initial and boundary conditions (|1.14p and p.!5p hold: 

- [ (v,ijH<f>)dt+ [ (v-Vv,^(p)dt 
Jo Jo 

-(vq, 0)^(0) + / (VdQVd,tpV(/))dt- (a,^V(j))dt, 
Jo Jo 

where a is defined in (jl.4p . and 

- / (d,ip t (f))dt+ / (v -X7d,ip(f))dt - / (nd,ip(/))dt + ^ / {Ad,%l)<p)dt 
Jo Jo Jo ^1 7o 



1 



T 



= -(do, 0)^(0) -—/ {Ad-f(d),^)dt. 
A i Jo 

With the help of the maximum principle under the assumption A2 = 0, in [32], the authors 
obtained the existence of weak solutions by applying a semi-Galerkin procedure (cf. [30] for the 
simplified liquid crystal system). For the more general case considered in the present paper, we 
no longer assume that A2 = 0. Consequently, the kinetic transport includes the stretching effect 
that leads to the loss of maximum principle for d. In order to ensure that the extra stress term 
V • a is well defined in the weak formulation (cf . Definition I4.ip , the regularity 

d G L°°(0,T; L°°) 

turns out to be essential (we refer to |44] for the discussions on the rod-like molecule liquid 
crystal model, which is a special case of the general system (|l.lip - (|1.15p ). In the subsequent 
analysis, we have to confine ourselves to the periodic boundary conditions, which helps us to 
avoid extra difficulties involving boundary terms when performing integration by parts in the 
derivation of higher-order energy inequalities. 

Finally, we remark that existence of global weak solutions to simplified liquid crystal systems 
in Remark 13.21 has been obtained in [3] with a suitable set of boundary conditions (i.e., homo- 
geneous Dirichlet boundary condition for v together with the homogeneous Neumann boundary 
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condition for d) . Their argument is based on an appropriate choice of test functions that leads 
to a suitable weak formulation of the system and thus overcomes difficulties from the stretching 
effect. Quite recently, the existence of global weak solutions with energy bounds to the general 
Ericksen-Leslie system (|l.lip - (|1.13|) has been proved in [3] by extending the argument in [3]. 
Moreover, in [4], under Parodi's relation, the authors prove the local existence/uniqueness of 
classical solutions to the general Ericksen-Leslie system (jl.llj) — (|1.13j) and establish a Beale- 
Kato-Majda type blow-up criterion. 



4.1 Galerkin approximation 

Under periodic settings, one can define a mapping S associated with the Stokes problem: Su = 
—Au for u G D(S) = {u £ H, Su G H} = Br? nil. The operator S can be seen as an unbounded 
positive linear self-adjoint operator on H. If D(S) is endowed with the norm induced by Hp 1 , 
then S becomes an isomorphism from D(S) onto H. 

Let {fa}^2 =1 with \\fa\\ = 1 be the eigenvectors of the Stokes operator S in the periodic case 
with zero mean, 

— Afii + VPj = Ki4>i, V • fa = in Q, / fa(x) dx = 0, 

JQ 

where P{ G L? and < k\ < ft 2 < ••• ar e eigenvalues. The eigenvectors fa are smooth and the 
sequence {fa} ( *L 1 forms an orthogonal basis of H (cf. |45j). Let 

P m : H ->• H m = span{(f>i, ■ ■ ■ , (f> m }, m £ N. 

We consider the following (variational) approximate problem: 

(d t v m ,u m ) + (v m ■ Vv m 

= (Vd m Vd m ,Vu m ) - (a m , Vu m ), V u m G H m , (4.3) 
A 2 1 

N m + —A m d m = -—Ad m - /(dm), (4.4) 
Ai Ai 

v m (x,0) = P m v (x), d m (x,0) = do(x), (4.5) 

v m (x + ei,t) = v m (x,t), d m (x + ei,t) = d m (x,t), (4.6) 

where 

n m = ^{Vv m - V T V m ), A m = ^{Vv m + V T i; m ), 

N m = d t d m + (v m ■ V)d m + tt m d m , 

a m = fii (d m Amd m )dm ®d m + ^N m ®d m + fj, 3 d m N m + ^iA m 
-\-\i^A m dm ® d m + n§d m ® A m d m - 

We can prove local well-posedness of the approximate problem (|4,3p - (|4,6|) by a similar semi- 
Galerkin procedure like [3] (see also [30 ] I32 | ). Smoothness of the approximate solutions in 
the interior of Qt = (0, To) x Q follows from the regularity theory for parabolic equations 
and a bootstrap argument (cf. |22,30j). The uniqueness of smooth solutions can be proved by 
performing energy estimates on the difference of two different solutions and using Gronwall's 
inequality. Since the proof is standard, we omit the details here. 

Proposition 4.1. Suppose that vq G V , do G Hp. For any m > 0, there is a Tq > depending 
on vq, do and m such that the approximate problem (|4.3j) — (j4.6|) admits a unique weak solution 
(v m ,dm) satisfying 

v m GL oo (0,r ;y)nL 2 (0,T ;H2), 
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d m €L 00 (0,r ;Hj)n£ 2 (0,To;Hj). 
Furthermore, (v m ,d m ) is smooth in the interior of Qt - 

4.2 Uniform a priori estimates 

In order to prove global existence of solutions to the problem fjl.Xlj) — f)1.15|> . we need some 
uniform estimates that are independent of the approximate parameter m and the local existence 
time To. These uniform estimates enable us (i) to pass to the limit as m — > oo to obtain a 
weak solution to the system (jl.lip - (|1.15p in proper spaces; (ii) to extend the local solution to a 
global one on [0, +00). One advantage of the above mentioned semi-Galerkin scheme is that the 
approximate solutions satisfy the same basic energy law and higher-order differential inequalities 
as the smooth solutions to the system (|1.11|) - (|1.15|) . For the sake of simplicity, the following 
calculations are carried out formally for smooth solutions. They can be justified by using the 
approximate solutions to (|4.3|) - (|4.6|) and then passing to the limit. 

The basic energy law plays an important role in the derivation of uniform estimates on 
L 2 x H 1 -norm of (v, d). According to the discussions in Section [21 we consider the following two 
cases, in which the basic energy law holds: 

• Case I (with Parodi's relation): Suppose A 2 / 0, ([LIB]) . (fTTTD . (|2~3j) -(l2il: 

• Case II (without Parodi's relation): Suppose A 2 / 0, (fTTBD . ([23|h([23D and (|2Tll . 
First, we consider Case I. It follows from Lemma [2JJ that 



4-£(t)<-[ Vi\d T Ad\ 2 dx-^\\Vv\\ 2 + ^-\\Ad- f(d)\\ 2 , Vt>0. 
at Jq 2 Ai 

This easily implies the following uniform estimates 

\\v(t,-)\\<C, |KV)[| H i <C, Vt>0, (4.7) 

J q °° (J Vi\d T Ad\ 2 dx + tj\\Vv\\ 2 -y\\Ad-f(d)\\ 2 ^Jdt<C, (4.8) 

where the constant C > depends only on ||uo|| and || cZo II hi - 

As we have mentioned before, the regularity d G L°°(0, T;L°°) is crucial to ensure that the 
extra stress term V • a can be suitably defined in the weak formulation. Due to the lack of 
maximum principle for d, an alternative way is to prove higher-order (spatial) regularity of d, 
e.g., in L°°(0,T;H 2 ) and use the Sobolev embedding H 2 <^-» L°° (n = 3). For this purpose, we 
derive a new type of higher-order energy inequality, which turns out to be useful in the study 
of global existence of solutions as well as the long-time behavior (cf. [30 , 32 , 44ll47j for simplified 
liquid crystal systems). 

Lemma 4.1. Set 

A(t) = \\Vv(t)\\ 2 + \\Ad(t)-f(d(t))f- (4-9) 

Let (J, be an arbitrary positive constant. We suppose that fj,^ > fj,. For n = 3, under the 
assumption of Case I, the following inequality holds: 

jA(t) < -{^-C x ^A{t))\\Avf 



+ \2X 1 + I l|V(Ad " 7)11 + C3A{t) ' (410) 



14 



where 

A{t) = A(t) + 1, 

Ci (i = 1,2,3,) are constants depending on Q, f, \\vq\\, H^oIIh 1 ; Ai, ^2, Hi (i = 1,2,3,5,6) and 
It- 
Proof. Without loss of generality, we assume that fi = 1. The argument is valid for arbitrary 
but fixed fi > 0. 

Using (|l.lip - (|1.13p and integration by parts, due to the periodic boundary conditions, we 
obtain that (see Section 17,31 for detailed computations) 

1 d f 

-— A(t) + m / {d k d p ViA kp fdx + y||At>|| 2 

J Q — 

Hhs + Hs) I IdjV^dx - h\V(Ad - f)\\ 2 
Jq a i 

= -fii / A kp V i(d k d p )didjV lAijdx - m / A kp d k d p V ' z (d;dj)V \A^dx 
JQ JQ 

— (M5 + Me) / ^ idjd k A ki V iA i: jdx - + n&) \ djV 'id k A ki V 'lA^dx 

JQ JQ 

- / Vj(Adi- f i )n ij Vid j dx+ / (Adi- f^VifLijVidjdx 

JQ JQ 

+2X 2 / N i ViA ij Vid j dx + \ 2 (N,AAd) 
Jq 

-^r~ I \^i{Aijdj)\ 2 dx + (Av,v -Vv) + -!- / f'(d)\Ad- f\ 2 dx 

-(Ad-/,/'(d)(Od-^4d)) +2 y VjiAdi- ftVwVtdidx 
-(Ad-f,vVf) 

/1 + ... + /14. (4.11) 



A 



In what follows, we estimate the right-hand side of (|4.1ip term by term. 

h<^ [ (dkdpV^pfdx + cwdWl^wvvWlsWvdWle. 

4 J Q 

By the estimate (|4.7p . we infer from the Agmon's inequality that 

WL»<C(l + ||Ad||*). (4.12) 

Then from (|4.7p . (|4.12p and the Gagliardo-Nirenberg inequality, we obtain 

[|Vu[| L 3 < \\Vv \\^\\Av ||Vu|| L 4 < [|Vi;||*[|Au[|i, (4.13) 

||Vd|| L 6<C(||Ad|| + l), (4.14) 
||Ad|| < ||Ad-/(d)|| + ||/(d)|| < ||Ad-/(d)||+C, (4.15) 

||VAd|| < ||V(Ad-/(d)|| + ||V/(d)|| 

< ||V(Ad-/(d)|| + ||/'(d)||L~ l|Vd|| 

< ||V(Ad-/(d)|| + C(l+||d||jU) 
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As a result, it holds 



< ||V(Ad-/(d)||+C(l + ||Ad||) 

< ||V(Ad - f(d)\\ + C(l + ||VAd||3||Vd||3 + ||Vd||) 

< ||V(Ad-/(d)|| + i||VAd||+C. 



\l~\\Vv\\l 3 \\Vdf L6 

< C||V<,||||A v ||(||Ad-/|| 3 + l) 

< Ll + 4 1| Ad - /|| 2 ) || A^|| 2 + C„J* ||V W || 2 (1 + || Ad - /|| 4 ) 

< 4^ll Au f + ^4*I|V«|| 2 

+C/i7l|Vvf (||VAd[|3||Vd[|3 + ||Vd|| +C 1 ) 4 

< ^A\\Av\\ 2 + C^\\Vvf + C M 4 ^||V«|| 2 (||V(Ad - /)|| 2 + 1) 

< 4A\\Av\\ 2 + Cfi4*A\\V(Ad - /)|| 2 + CA, 
which implies that 

h < ^ f {didjViAitfdx + 4A\\Av\\ 2 

4 JQ 

+Cfi~^ A\\V(Ad - f)\\ 2 + CA. 
For I2, using integration by parts, we obtain 

h = Mi / ViA hp d k d p Vi(didj)Aijdx + m / A kp Vi{d k dp)Vi{didj)Aijdx 
JQ JQ 

/ A kp dkd p (djAdi + 2VidiV[dj + diAdj)A^dx 

JQ 

< Y f (4d P V^ fep ) 2 dx + C||d|| 2 j oo||V W || 2 3 ||Vd|| 2 6 



where 



+C||V V || 2 4 ||Ad||||d||i co , 



C\\Vv\\l4Ad\\\\d\\l a 



< C||Vu||2||Au||2(||Ad-/||2 +1) 

< M 4 U||Aw|| 2 + Cfi 4 -i\\Vvf(l + || Ad - /|| 4 ). 

Thus, the right-hand side of (|4.20p can be estimated exactly as (|4.17p . Therefore, 

h < y / (4d P V^ fep ) 2 dx + / u|^l|Av|| 2 

+C/J 4 M||V(Ad-/)|| 2 + CA 
Concerning ^3 and I4, we deduce from (|4.16p that (using again (|4.17p and fj,^ > 1) 

I3 + /4 = -(^5 + ^6) / V 'idjd k A ki V lAijdx 

JQ 



IQ 
16 



-(^5 + /^) / djV id k A ki V lAijdx 
JQ 

< C||A W ||||Vt;|| L 3||Vd|| L6 ||d|| L ^ 

< /i45||At;f + CM4-5||Vw||£a||Vd||?,«||d||f l00 

< vMM? + C\\Vv\\l4Vd\\l4d\\l.o 

< 2fijA\\Av\\ 2 + C^A\\V(Ad- f)f + CA. (4.22) 



Next, for 1$, Iq, we have 



h < q|v(Ad-/)||||v w || L3 ||Vd|| L 6 

< C||V(Ad-/)[|||Vt;||3||At;||3(||Ad-/[| + l) 

< vi\\Vv\\\\Av\\+C»~*A\\V(Ad-f)f 

< fil\\Avf + C^A\\V(Ad-f)f + C\\Vv\\ 2 , (4.23) 



and 



h < ||vn||||Ad-/||||vd|| L ~ 

< C||A w ||||Ad-/||(||V(Ad-/)||f + 1) 

< 4l|A^|| 2 + C^ § ||Ad-/|| 2 ||V(Ad-/)|| 2 + C||Ad-/|| 2 

< ni\\Avf + Cn^A\MAd-f)\\ 2 + CA. (4.24) 
Using integration by parts and (|1.13p . we get 

h + h + h 

= 2A 2 f NiVtAijVidjdx + X 2 (N,AAd) - [ \V ^d^ 2 dx 

JQ A l JQ 

= A 2 / NiViAijVidjdx- A 2 / V 'M^V \djdx 



Ai 

(A2) 2 
Ai 



(A ^- / lAijVtdtfdx-^- [ VtAijdjAuVidkdx 
*i JQ A i JQ 

— [ iViAijdjfdx 

^ Jq 

[ (Adi- f i )V l A ij Vid j dx+^- [ Vi(Adi- f^AijVidjdx 
M Jq M Jq 

Al Jq Ai J q 

K 1 +K 2 + K 3 + K A . (4.25) 



Similar to (ET221) . (1031 and P~24"D . we have 



< C||At;||||Ad-/||||Vd||L- + C||V(Ad-/)||||V«|| t 3||Vd|| L fl 

< 4 ||A«|| 2 + C^A\\V(Ad - f)\\ 2 + CA, 

K 3 < fxjA\\ Avf + C^ A -'i\\Vv\\ 2 . 
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Furthermore, (|2.4p and (|2,6p indicate that 

K4 - (ji 5 + Me) / iV^-djfds < 0. (4.26) 

JO 



As a result, 



I 7 + h + h 

< hZA\\Av\\ 2 + (h 5 + hg) / {ViAijdjfdx 

Jq 

+C^A\\V{Ad- f)\\ 2 + CA. 
For terms I%q, In and /12, we have 

ho < ||u|| L 4||Vu|| L 4||Au|| 

< C||«||? ||V«||^ [|Vv||i||Ai;||t \\Av\\ 

< ^\\Avf + n^\\w\\ 2 \\Av\\ 2 + C^-^WVvf 

< fi^A\\Av\\ 2 + CA, (4.27) 



111 < C(\\d\\le + l)\\Ad- f\\l 3 

< c(||Ad-/||||V(Ad-/)|| + ||Ad-.f|| 2 



and 



< -±-\\V(Ad-f)\\ 2 + C\\Ad-ff, (4.28) 



h2 = -(Ad-f,f{d)(nd-^Ad 

< C(||d||= 8 + l)||d|| Le ||Ad-/|| La [|Vt;[| L 6 

< c(||V(Ad- /)|| + ||Ad- /||)||Vu||*||Ai;||* 

< 4ll Aw ll 2 + ^^||V(Ad-/)|| 2 + CA (4.29) 
The estimate for /13 is exactly the same as (|4.23|) such that 

J 13 < ^||A«|| 2 + C^A\\V(Ad - f)\\ 2 + C||V W || 2 . (4.30) 
Finally, for /14, we see that 

I 14 < C||Ad-/|| L 3H| L 6||V/|| 

< ciiAd-ziiLaiiv^iKi + Hdii 2 ^)!^^! 

< C(l + ||Ad - /||)(||V(Ad - /)|| + II Ad - f\\)\\Vv\\ 

< + (-4^ + vt 4 ) HV(Ad - /)|| 2 + C(l + H ^)A. 

Putting all the above estimates together, we arrive at the higher-order differential inequality 
(|4.10p . The proof is complete. □ □ 
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di A{t) + \2 -CiWA(t))\\Av\ 



Lemma 4.2. Under the assumption Case I, for any initial data (yQ,do) £ V x H 2 , if the 

viscosity \i± is properly large, i.e., 

(14 > nl(jj, i ,\i,\2,vo,do,tJ L ), i = 1,2,3,5,6, 

we have 

A(t)<C, Vt>0. (4.31) 

The uniform bound C is a constant depending only on f , Q, ||t>o||y, H^oIIh 2 ; I 1 ' 8 ' ^' s > V" 

Proof. It follows from (|4.10p that 

d 
Tt 

+ ("^ - C 2H l A(t)^ ||V(Ad - /)|| 2 
< C 3 A(t). (4.32) 
Meanwhile, by (|4.8p . we have 

/ A(r)dT < / „4(r)dr + 1 < M, Vi > 0, (4.33) 

where M is a positive constant depending only on /^s (except ^4), A^s, ||«o||> II^oIIh 1 - Now we 
choose H4 large enough such that 

fij > 2d (A(0) + 4M + C 3 M) + 4A 2 Cf (.1(0) + 4M + C 3 M) 2 + 1. (4.34) 

Applying a similar argument in [32, Theorem 4.3] (cf. also [301IIE]), we can see that A(t) is 
uniformly bounded for all £ > and satisfies 

^ - C lfJ ,^A(t) > 0, -1- - C 2 ^A(t) > 0, Vt>0. (4.35) 
2 — 2Ai 

The proof is complete. □ □ 

Next, we briefly discuss Case II. 

Corollary 4.1. For n = 3, under the assumption Case II, the inequality (|4.10p still holds. 

Proof. If Parodi's relation (|l,17p does not hold, i.e., A2 + (fJ-2 + ^3) / 0, then in the derivation 
of -^A(t) (see Section [T7S|) . the first term on the right-hand side of (|7.19|) does not cancel with 
the first term on the right-hand side of (|7.16p . Consequently, there is one extra term: 



(A 2 + /u 2 + ^3) / dj^AA^dx. 
Jo 



Besides, since we no longer have (|2.6p in Case II, we have to re-investigate the left-hand side 
of (|4.26p . Using the d equation f|l . 13j) and integration by parts, we get 



(A 2 + H2 + M3) / djNiAA^dx 
Jo 



A2 + fJ-2 + ^3 



/ djVi{Adi - fi)ViAijdx 
Jq 

+ A 2 + /i 2 + ^ 3 f Vldj{Adi _ fi)VlAjdx 
M Jo 
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+ 



+ 



+ 



A 2 (A 2 + /u 2 +/U3) 





Ai 


A 2 (A 2 


+ + yU3) 




Ai 


A 2 (A 2 


+ ^2 + ^3) 



|djVi^4ij| 2 drc 



Q 



VidjAikdkViAijdx 



Ai 



Q 



djAikVidkViAijdx. 



(4.36) 



We now estimate the right-hand side of (|4.36p . For the first term, we have 

A 2 + /u 2 + M3 



Ai 



/ djVi(Adi - fi)ViAijdx 
Jo 



< C||d|| L oo||V(Ad-/)||||Au| 

< C(||Ad-/||s+l)||V(Ad- 

< ^(1 + \\Ad - f\\)\\Av\\ 2 - 



Av\ 



C 



r ||V(Ad-/)|| s 



M4 4 

|2 i ri.rt IIV7/A J *M|2 



< M4 5 ^||A V r + C/i 4 4 ||V(Ad-/)||^ 



(4.37) 



The second term can be estimated as (|4.24p . while the fourth and fifth terms are similar to 
(|4.22p . Finally, concerning the third term and the two terms on the left-hand side of ()4.26p . we 
infer from ([1716]) and (f2"7IT]) that 



A 2 (A 2 + M2 + M3) (A2 



Ai 



Ai 



(^5 + A*6) 



■ [Ai(//5 + Me) - A 2 (^2 + M3)] 



K "IT 
1 



-(A 2 - a*2 - M3) 2 - A 2 (m2 + M3) 



^[(A 2 ) 2 + (/i 2 + ^ 3 ) 2 ] <0, 



which yields 



A 2 (A 2 + ^ 2 + ^ 3 ) (A 2 ) 5 



(A*5 + Me) 



djVi^rda; < 0. 



Ai A x 

Combining the other estimates in the proof of Lemma [4.1] we obtain the inequality (|4.10p under 
assumption Case II. □ □ 

Corollary 4.2. Under the assumption Case II, for any initial data (vo,do) 6 V x H 2 , if the 

viscosity \i± is properly large, i.e., 

M4 > fJ-lifH, Ai,A 2 ,uo,do,//), « = 1,2,3,5,6, 

we Ziaue «4(i) < C for t > imi/t C feeing a constant depending only on f, Q, \\vq\\v, ||^o||h 2 > 
H 1 s, X's and [i. 

4.3 Global existence and uniqueness 

In both Case I and Case II, the uniform estimates we have obtained in Section 4.2 are in- 
dependent of the approximation parameter m and time t. This indicates that for both cases, 
(v m ,d m ) is a global solution to the approximate problem f|4.3j) — (|4=.5j) : 



v m G L°°(0,+oo;y)nL 2 oc (0,+oo;H 2 ), 
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d m G L°°(0,+oo;H2)nLL(0,+oo;H3), 

which further implies that 

d t v m G Lf oc (0, +00; Lj), G Lf oc (0, +00; Hj). 

The uniform estimates enable us to pass to the limit for (y m ,d m ) as m — > 00. By a similar 
argument to [30, ,33], we can show that there exist a pair of limit functions (v,d) satisfying 

v G L°°(0, 00; y) n Lf oc (0, +00; Hj), (4.38) 
d G L°°(0, +00; Hp) n L 2 OC (0, +00; Hp), (4.39) 

such that (v,d) is a weak solution of the system (|l.lip - (|1.15|) . A bootstrap argument based on 
Serrin's result |42] and Sobolev embedding theorems leads to the existence of classical solutions. 
The uniqueness of solutions to the problem (|1.11|) - (|1.15|) with regularity (|4.38p - (|4.39p can be 
proved as in [451 Lemma 2.2]. 
In summary, we have 

Theorem 4.1 (Global well-posedness). Let n = 3. We assume that either the conditions in 
Case I or in Case II are satisfied. For any (vo,do) G V x Hp, under the large viscosity 
assumption 

fj,4 > n1(jii,\i,\2,vo,do,fi), i = 1,2,3,5,6, 
the problem fjl.Xlj) — f|l. 15j) admits a unique global solution that satisfies (|4.38p - ()4.39|) . 

Besides, we have the following continuous dependence on the initial data: 

Lemma 4.3. Suppose that the assumptions in Theorem \4-l\ are satisfied. {yi,di) {i = 1,2) are 
global solutions to the problem (|1 . 11 (j 1. 15 j) corresponding to initial data (foi,c?oi) £ V x Hp 
(i = 1, 2). Then for any t G [0, T], we have 

\\(vi-V2)(t)\\ 2 + \\(d 1 -d 2 )(t)\\ 2 lll 

+ f (yl|V(«i - ^ 2 )(r)|| 2 + ||A(di - d 2 )(r)f) dr 
< 2e ct {\\v 01 - v 02 \\ 2 + \\doi - doallni). 
where C is a constant depending on ||«o*||V) II^OjIIh 2 ) ^' s ) but not on t. 
Remark 4.1. If in addition, we assume either 

(i) fix = 0, A 2 / 0, or (ii) m > 0, A 2 = 0, 

the same result holds true in 2D without the largeness assumption on /X4. In case (i), we note 
that the nonlinearity of the highest-order vanishes. In particular, this applies to the system 
(|3.14p - (|3.16p . which is a simplified version of the general Ericksen-Leslie model (cf. Jj4tf for the 
liquid crystal system with rod-like molecules and f^ [To\\48\l with general ellipsoid shape). On the 
other hand, in case (ii), one can apply the maximum principle for d to obtain its L°° -bound, 
which makes the proof much easier (cf. JEEj)- D 
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4.4 Long-time behavior: convergence to equilibrium 

Now we briefly discuss the long-time behavior of the global solution (y, d) obtained in Theorem 
14.11 First, we have the following decay property: 

Lemma 4.4. For the global solutions obtained in Theorem\4-l\ we have 



lim (K*)|| v + ||-Ad(t) + /(d(t))||) = 0. (4.40) 

t— >+oo 

Proof. We only consider Case I and the proof for Case II is similar. From the basic energy 
law (|2.7p . we see that A(t) G L^+oo). On the other hand, (|4TTnD together with (OT]) and 
(|4.35|) implies that 4iA(t) < C. As a consequence, 

lim A(t) = 0. 

i->+oo 

The proof is complete. □ □ 

It easily follows from Lemma 14.41 that 
Proposition 4.2. Suppose that the assumptions in Theorem \4-l\ are satisfied. The u -limit set 



of (vo,do) £ V x Hp denoted by cj(i>o,do) is a non-empty bounded connected subset in V x H 2 , 
which is also compact in L 2 x EL. Besides, we have 

oo(vq, do) G S := {(0, d) : —Ad + f(d) = 0, in Q, d(x + &i) = d{x) on dQ}. 

Therefore, all asymptotic limit points of the system ([l.lip — (|1.13p satisfy the following reduced 
stationary problem 

Voo = 0, (4.41) 
VPoo+vf^^ = -Vcioo-Adoo, (4.42) 



2 

-Adoo + /(doo) = 0, (4.43) 
d^x) = dooix + a), xedQ, (4.44) 



where in (|4.42p . we used the well-known fact that (cf. [30] ) 

V • (Vdoo Vdoo) = V f l W o °°l ) + Vdoo • Adoo. 



(|4.42p is a constraint equation for doc. Since T'[d) = f(d), if doo is a solution to (|4.43p . then 
(|4.42p is automatically satisfied because all gradients can be absorbed into the pressure. 

We have already proved that the velocity field v decays to zero in V as t /* +00 (cf. Lemma 
I4.4p . On the other hand, we can only conclude sequential convergence for d from compactness 
of the trajectory: for any unbounded sequence {tj}, there exist a subsequence {t n } /* +00 such 
that 

lim ||d(t n )-d 00 || H i =0, (4.45) 

where d^ satisfies ()4.43p - (|4.44p . The convergence of d for the whole time sequence is non-trivial 
because in general we cannot expect the uniqueness of critical points of E(d). In the present case, 
under the periodic boundary conditions, one may see that the dimension of the set of stationary 
solutions is at least n. This is because a shift in each variable may give another steady state. The 
convergence of d to a single equilibrium can be achieved by using the well-known Lojasiewicz- 
Simon approach (cf. L. Simon [43]). We refer to [19] and the references therein for various 
generalizations and applications. To this end, we introduce a suitable Lojasiewicz-Simon type 
inequality in the periodic setting (cf. e.g., [I9~]). 
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Lemma 4.5 (Lojasiewicz-Simon inequality). Let ip be a critical point of the functional 

E(d) = -\\Vd\\ 2 + / F{d)dx. 

Then there exist constants 9 € (0, ^) and /3 > depending on ip such that for any d € 
satisfying \\d — V'IIh 1 < fl, it holds 

||- Ad + /(d)|| (H iy > \E(d)-EW\ 1 - e , (4.46) 

where (H*)' is the dual space of Hp. 

Then we have the following convergence result: 

Theorem 4.2 (Convergence to equilibrium). Under the assumptions of Theorem \4- 1\ the global 
solution (v, d) has the following property: 

IKt)||v + IKt)-doo||H? <C(l + t)~(&V, Vt>0, (4.47) 

where doo is a solution to (|4.43p ~ (|4.44p . C is a constant depending on vq, do, f, Q, fi^s, \[s, 
doo and the constant 6 6 (0, ^) depends on d^ (called Lojasiewicz exponent, cf. Lemma \4.5ty . 

Based on Lemma 14.51 the basic energy law deduced in Section 2 (cf. (|2.7p or (|2,12p ) and 
the higher-order energy inequality (cf. Lemma 14.11 or Corollary I4.ip . we can prove Theorem 14.21 
following the procedure in [481 Section 3.2, 3.3] with minor modifications. In order not to make 
the paper too lengthy, we leave the details to interested readers. 



5 Well-posedness and Nonlinear Stability under Parodi's Rela- 
tion 

The results obtained in Section U] indicate that for both Case I (with Parodi's relation) and 
Case II (without Parodi's relation), global well-posedness of the Ericksen-Leslie system can be 
obtained provided that the viscosity fi^ is properly large. Recall the Navier-Stokes equations in 
3D (with periodic boundary conditions and vq G H), we can easily derive 

J t \\Vv\\ 2 + (^4-/i45||V V || 2 )||At;|| 2 < C^^IIV^H 2 , 

which implies that the large viscosity assumption is equivalent to small initial data assumption on 
v in H 1 -norm. However, this is not the case for the Ericksen-Leslie system f|1.15|) due to its 

much more complicated structure (cf. (|4.34p ). Actually, we do not have the large viscosity /small 
initial data alternative relation even for those simplified liquid crystal systems [30 1144] . 

In this section, we show that Parodi's relation (|1.17p plays an important role in the well- 
posedness and stability of the system (|l,lip - (|1.15p . if no additional requirement is imposed on 
the viscosity [14. In particular, under those assumptions in Case I, we are able to prove a 
suitable higher-order energy inequality that yields the local well-posedness and furthermore, the 
global existence result provided that the initial velocity vq is near zero and the initial director 
do is close to a local minimizer d* of the elastic energy 



E(d) = l\\Vd\\ 2 + [ F{d)dx. (5.1) 

2 Jq 



Besides, we are able to show the Lyapunov stability of local energy minimizers of E(d). This 
implies that Parodi's relation (|1.17|) serves as a sufficient condition for nonlinear stability of the 
Ericksen-Leslie system (jl.lip - (|1.15p from the mathematical point of view. 
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5.1 Higher-order energy inequality and local well-posedness 

Lemma 5.1. Letn = 3. Suppose that the conditions in Case I are satisfied. Then the following 
higher-order energy inequality holds: 



j t A(t) + ^ J (d^V^fdx + ^Avf-^-llViAd-fn 2 



>Q 

< C*(A 6 (t)+A{t)), (5.2) 



where C* is a constant that only depends on (j! s, X's, \\vq\\ and \\do\ 



H 1 



Proof. First, from the basic energy law (|2.T|) we still have the uniform estimates on ||f(t)|| and 
II^WIIh 1 Q4.7p ). Moreover, estimates f|4. 12|) — (|4. 16[) are still valid. Next, we re-estimate the 
terms I\ , . . . , /14 on the right-hand side of (|4.1ip . 

h < ^ / (dkdpVtA^fdx + C\\d\\iao II Vv\\ls\\Vd\\le. 

< ^ / (dkdpV^fdx + CiWAd-ff + l^VvWWAvW 

< a f (d^V^fdx + ^-WAvW 2 

+C(\\Vv\\ 2 \\Ad- ft + \\Vvf) 

< ^7 f {d k d p ViA kp ) 2 dx + ^\\Avf + CA A + CA, (5.3) 

4 Jq ol 

h < ^ I (d k dpViA kp ) 2 dx + C||d||L»||V«||£a||Vd|£ 8 

+q|d||f,^V< 4 ||Ad|| 

< ^ / (4d P V^ fcp ) 2 (ix + ^||Ar;|| 2 + C7^ 4 + C7^ 

+C'(||Ad-/||3 + l)||V«||3||At;||3 

< !j f (d k d p ViA kp ) 2 dx + ^\\Avf + CA 6 + CA, 



h + h < C||A V ||||V«|| L 3||Vd|| L 6||d|| L - 

< gllA^f + CUVHlLllVdlllelldll^ 

< ^\\Av\\ 2 + CA 4 + CA, (5.4) 



h < C||V(Ad-/)||[|V«[|3||At;[|5([|Ad-/[| + l) 

< gl|A^|| 2 - ^l|V(Ad - /)|| 2 + C||Vt,|| 2 (||Ad - /|| 4 + 1) 

< ^\\Av\\ 2 -^-\\V(Ad-f)\\ 2 + CA 3 + CA, (5.5) 



h < q|A W ||||Ad-/||(||(V(Ad-/)||4+l) 
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< g||A W || 2 -^-||V(Ad-/)|| 2 + ^ 4 + CA (5.6) 

For the terms K\, ...,^4 in (|4.25p . we still have ()4.26p . By a similar argument to f)5.5|) f)5.6[) we 

get 

K, < !^\\Avf-^-\\V(Ad-f)f + CA' + CA, 
K 2 < ^\\Av\\ 2 -^-\\V(Ad-f)f + CA" + CA. 
Using integration by parts, we can see that 

K 3 = -^M- f ViAijdjAijVtdjdx 
A i Jq 

= ( \A ij V l d j \ 2 dx-K 3 + ^- f AijdjAijAd 

M Jq M Jq 

which together with similar estimates in fj5.4[) yields that 

K 3 = ^£ / \A ij V l d j \ 2 dx + / AijdjAijAdjdx 

Ai Jq Ai Jq 



Hence, 



< ^\\Av\\ 2 + CA* + CA. 



I 7 + h + h 



< ^\\Av\\ 2 - ^||V(Ad - f)\\ 2 + 0*6 + Me) / |V,^,^| 2 cfa 
+CA 4 + CA 

The remaining terms can be estimated in a straightforward way. 

iio < \{Av,v -Vv)\ <C\\Av\\ 7 i\\Vv\\ 
< g||A,|| 2 + C7||V,|| 8 , 



In < C(||d|| 2 6 + l)||Ad-/|| 2 3 
< --L\\v(Ad-f)\\ 2 + CA, 



I12 < C\\f'(d)d\\\\Ad-f\\- L 4Vv\\- L3 

< C||V«||3||A«||^||V(Ad-/)|| + ||Ad-/|| 

< ^\\Av\\ 2 -^-\\V(Ad-f)\\ 2 + CA. 



The estimate of I\ 3 is similar to (J53j): 



hs < - g^-||V(Ad - f)\\ 2 + CA 3 + CA. 
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For the last term Z14, we have 



I U < C||Ad-/|| L 3|H| L6 ||V/|| 

< C(l+||Ad-/||)(||V(Ad-/)|| + ||Ad-/||)||V«|| 

< -^-\\V(Ad-f)\\ 2 + CA 2 + CA. 

0A1 

Collecting all the estimates above, we can conclude the higher-order differential inequality (|5.2p . 
The proof is complete. □ □ 

The following local well-posedness result is a direct consequence of the higher-order energy 
inequality (|5.2|) : 

Theorem 5.1 (Local well-posedness). Let n = 3. Suppose that the conditions in Case I are 
satisfied. For any (vo,do) G V x H 2 , there exists a T* > such that the problem f|l. 11 j) — t|l. 15j) 
admits a unique local solution satisfying 

v e L°°(0,T*-V)r\L 2 (0,T*;Ul), d£ L°° (0, T*;H 2 ) n L 2 (0, T* ; H^) . 

Remark 5.1. Unfortunately, we are not able to prove a corresponding local well-posedness result 
under the assumptions in Case II where Parodi 's relation (|1.17p is not satisfied. In this case 
the higher-order energy inequality (|5.2p is not available any longer. One obvious difficulty is that 
we lose control of some higher-order nonlinearities that will vanish due to specific cancellations 
under Parodi 's relation (see, e.g., (|4.37|) ). □ 

5.2 Near local minimizers: well-posedness and nonlinear stability 

Based on Lemma 15, l\ one can easily deduce the following property: 

Proposition 5.1. Suppose that the assumptions in Case I are satisfied. For any (vo,do) £ 

v x n 2 p , if 

||Vv|| 2 (0) + II Ad - /(d)|| 2 (0) < R, (5.7) 

where R > is a constant, there exists a positive constant Sq depending on fjfs, X's, \\vq\\, 
II^oIIh 1 ? f> Q an d R> such that the following property holds: for the (unique) local solution (v,d) 
of the system (| 1 . 1 1 j) — fj 1 . 1 5 j) which exists on [0, T*\, if 

£(t)>£(0)-s , Vie[0,T*], 

then the local solution (v,d) can be extended beyond T* . 

Proof. We consider the following initial value problem of an ordinary differential equation: 

^y(t) = C*(Y(tf + Y(t)), Y(0) = R>A(0). (5.8) 
at 

We denote by / = [0,T max ) the maximal existence interval of Y(t) such that 

lim Y(t) = +00. 

It follows from the comparison principle that for any t £ I, < A(t) < Y(t). Consequently, A(t) 
exists on /. ^Ve note that T max is determined by y(0) — R and such that T max — T max (R, C*) 
is increasing as R decreases. Taking 

3 

U) = -rT max {R, C*) > 0, 
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then we have 

< A(t) < Y(t) < K, Vt € [0, t ], (5.9) 

where K is a constant that only depends on R,C*,to- This fact combined with the Galerkin 
approximate scheme in Section 4.1 leads to the local existence of a unique solution to the system 
(|l.lip - (|1.15p at least on [0,to]. (This indeed provides a proof of Theorem 15.11 ) 

The above argument suggests that the existing time T* > to- Now if £ (t) > £ (0) — Eq f° r a U 
t G [0, T*], we infer from Lemma 12. II that 

'^|V*;(t)| 2 - y \Ad(t) - f(d(t))\ 2 ) dxdt < e . 
Hence, there exists a t* G [T* — |,T*] such that 

||V^*)|| 2 + ||Ad(t,) - f(d(U))\\ 2 < max{ — 
Choosing £o > such that 

maxj — = R, (5.10) 

we have «4(t*) < -R. Taking t* as the initial time and Y{t*) = R in (|5.8p . we infer from 
the above argument that Y(t) (and thus A(t)) is uniformly bounded at least on [0,t* + to] D 
[0,T* + |t ]. Thus, we can extend the local solution (v,d) from [0,T*] to [0, T* + |t ]. The 
proof is complete. □ □ 

Remark 5.2. Proposition 15. 1\ implies that, for the local solution (v,d) of (|l.lip - (|1.15p . if the 
total energy £(t) does not drop too much on its existence interval [0,T*] ; then it can be extended 
beyond T* . We note that stronger results have been obtained in [30, 32} for simplified liquid 
crystal systems. In those cases, global existence of weak solutions can be proved and the total 
energy £(t) is well-defined on [0, +oo). Then one can show the alternative relation: either there 
exists a T < +oo such that £{T) < £ (0) — Eq or the system admits a (unique) global strong 
solution. □ 

Remark 5.3. It is easy to verify that the above hypothesis on the changing rate of £(t) can be 
fulfilled, if the initial velocity vq is near zero and the initial molecule director do is close to an 
absolute minimizer of the elastic energy E(d) (for instance, a constant vector with unit length). 
We refer to }30[\32[ \4 7 ^ for the cases of the simplified liquid crystal system. The same result 
holds for our current general case if the same assumption is imposed. □ 

The assumption that the initial director do is close to an absolute energy minimizer can 
indeed be improved. Under Parodi's relation (|1.17p . we can show much stronger result that if 
vo is near zero and do is close to a local minimizer of E(d), then the total energy £ will never 
drop too much. Actually, we shall see below that the global solution will stay close to the given 
minimizer for all time (i.e., Lyaponov stability) and £(t) will converge to the same energy level 
of the local minimizer. This generalized result also applies to all those simplified Ericksen-Leslie 
systems considered in the literature [3T) l l3"2" | l4"4 4l47l l4"8] . 

Definition 5.1. d* € is called a local minimizer of E(d), if there exists a > 0, such that for 
any d £ satisfying \\d — d*\\ni < o~, it holds E(d) > E(d*). 

Remark 5.4. Since any minimizer of E(d) is also a critical point of E(d), it satisfies the 
Euler-Lagrange equation 

-Ad + f(d) = 0, xeQ, d(x) = d(x + ei), x e dQ. (5.11) 

From the elliptic regularity theory and bootstrap argument, one can easily see that if the solution 
d G then d is smooth. □ 
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Next, we state the main result of this section: 

Theorem 5.2. Suppose that n = 3 and the conditions in Case I are satisfied. Let d* G Hp be 
a local minimizer of E{d). There exist positive constants o~i,o~2, which may depend on X^s, fi^s, 
Q, a and d* , such that for any initial data (vo, do) G V x H 2 satisfying 

\\v \\ n i < 1, ||d - d* || H 2 < 1 

and 

\\vq\\ < cri, || d - d*|| H i < (T 2 , 

we have 

(i) the problem () 1 . 1 1 p — fj 1 . 1 5 j) admits a unique global solution (v,d), 

(ii) (v,d) enjoys the same long-time behavior as in Theorem In addition, 

lim 8{t) = Efdoo) = E(d*). (5.12) 

t— >+oo 

Proof. Without loss of generality, we assume that the constant a in Definition 15 . 1 1 satisfies a < 1. 
Throughout the proof, Cj, i = 1,2, • • • denote generic constants depending only on fj^s, X^s, a 
and d* . By our assumptions, we easily see that 

Mt)\\ + ||d(t)|| H i <Cx, Vt>0, (5.13) 
A(0) = \\Vv \\ 2 + \\Ad - f(d )\\ 2 < C 2 . (5.14) 

Recalling the proof of Proposition ^. 11 we take R = C2 for our current case. The constant C* 
in (|5.8p can be determined by C\ and fx's, X's (cf. Lemma IBTT]) . Then we set to = \T max (C2, C*) 
and take T* = to- Finally, the critical constant £q is given by (|5.10p . It follows from (|5,9p that 
A(t) is uniformly bounded on [0, to], which implies 

Ki)||y + ||d(t)|| H 2 < C 3 , Vi G [0,to]. (5.15) 

Next, we extend the local solution to [0, +00) by using the Lojasiewicz-Simon approach. 
Since the minimizer d* is a critical point of E(d), we take ip = d* in the Lojasiewicz-Simon 
inequality (cf. Lemma H~5j) . then the constants /3 > 0, 6 G (0, 5) are determined by d* and (I4.46|) 
holds. 

The proof consists of several steps. 

Step 1. In order to apply Proposition 15. II with T* = to, it suffices to show that 

£ (t) - £ (0) > -e„, Vt€[0,to]. (5.16) 
Using (|5.13p and the Sobolev embedding theorems, we have 

\E(d ) - E(d*)\ <C 4 ||do-(f || H i, 

which implies that 

£(t)-£(0) = \\\v{i)\\ 2 -\\\vo\\ 2 + E{d(t))-E{do) 

> ~\\v \\ 2 + E(d(t)) - E{d*) + E{d*) - E{do) 

> -\\\vof -C4do-d*\\ m +E{d{t))-E(d*). (5.17) 



Take 



at < min|e 2 ,l|, ct 2 < min l|. 
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Then by (|5.17p . it is easy to check (|5.16p will be satisfied provided that 

E{d(t)) - E(d* ) > 0, Vt€[0,io]. (5.18) 
By the definition of d* , it reduces to prove that 

\\d(t) -d*\\ H i < a, ViG[0,i ]. (5.19) 
Actually, we can prove a slightly stronger conclusion such that 

\\d(i)-d*\\ H i <lj:= -min{o-,/3}, ViG[0,t o ]- (5.20) 

Suppose 

1 

We use a contradiction argument. If (|5.20p is not true, then by the continuity of d that d G 
C([0, to]; H 1 ), there exists a minimal time To G (0,to] ; such that 

\\d(T ) - d*\\ n i = gj. 

Observe that 

£(t) = \\\v(t)\\ 2 + E(d(t)) > E(d*), Vt G [0,T ]. 

First, we consider the trivial case that for some T < Tq, £(T) = E(d*). Then we deduce 
from the definition of the local minimizer that for t > T, £ cannot drop and will remain E(d*). 
Thus, we infer from the basic energy law (|2.7p that the evolution will be stationary and the 
conclusion easily follows. 

In the following, we just assume £(t) > E(d*) for t G [0,To]. Applying Lemma [4.51 with 
ip = d*, we get 

~[£{t) - E(d*)} 6 
= -9[£{t) - E(d*)f-^ t £(t) 

e{^\\Vvf-^\\M-ff 

- C{\\v\\ 2 ^- ) + \\Ad- f\\) 

> C 5 (||V^|| + ||Ad-/||), VtG(0,T ). 

On the other hand, it follows from fjl. 13j) and (|5.15p that 

Udell < \\v-Vd\\ + \\nd\\+ ^\\Ad\\-h\Ad-f\\ 

< C6(lkl| L 6||Vd|| L 3 + ||Vw||||d||i,=c + ||Ad - /||) 

< C 7 (||Vt;|| + ||Ad-/||) J V*€[0,t ]. (5-21) 



Consequently, 



\\d(T ) - doWui < C 8 \\d(T )-dop\\d(T )-do\\l 



i 

2 



< c 9 (£° \\d t (t)\\dty < c w [£(o) -E(d*)}-2 

< Cio(^KI| 2 + C4||do-cf||Hi V 
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< CuONI* + \\do-d*\\* ll ). (5.22) 
Finally, choosing (also taking the previous assumptions into account) 

o 1 = mm{4.(4^)'»l}. ^^^{^(i^) 9 '!' 1 }' (5 - 23) 
we can deduce from (|5.22j) that 

\\d(T ) - d*\\ u i < \\d(To)-do\\ n i + \\do-d*\\ H i 

U OJ U) 

* 4 + 4 + 4 <W ' 

which leads to a contradiction with the definition of To. Thus, (|5.20p is true and so is (|5.18p . 
which implies that (|5. 16[) is satisfied. 

As in the proof of Proposition 15.11 there exists a i* E [^,to], such that A(t*) < R. Then 
we conclude that A(t) is uniformly bounded on [0, t* + to] D [0, ^] (with the same bound as on 
[0, to]). Here, we note the important fact that the bound of A(t) only depends on R,C*,to but 
not on the length of existence interval. 

Step 2. Now we take T* = |io- By the same argument as in Step 1, we can show that 

£ (t) - £ (0) > -£ , t€[0,T*]. 

Again, we obtain that A{t) is uniformly bounded on [0, T* + |io] (with the same bound as on 
[0, to])- By iteration, one can see that the local solution can be extended by a fixed length |io 
at each step and A(t) is uniformly bounded by a constant only depending on R, C*, to- 

Therefore, we can show that (v, d) is indeed a global solution. Moreover, the following 
uniform estimate holds 

Kt)|| H i + ||d(t)||H?<^ Vi>0, (5.24) 

where K depends on Ci, R,C*,to- The conclusion (i) is proved. 

Step 3. Based on the uniform estimate (|5.24p . a similar argument to Theorem 14.21 yields 
that there exists a d^ satisfying ()4.43p - (|4,44p . such that 

lim (|Kt)|| v + \\d(t) - cUIhO = 0, (5.25) 

t— >+oc 

with the convergence rate (|4.47p (We remark that in (|4.47p . the Lojasiewicz exponent 9 is 
determined by the limiting function doo, which is different from the one we have used in Step 
1). 

By repeating the argument in Step 1, we are able to show that 

\\d(t)-d* || H i <u, Vt>0. 
Then for t sufficiently large, we have 

Moo-^||hi < ||d^o — d(t)||ni + ||d(t) — tT[|iri 

< \j < mm{/3,a}. (5.26) 

Applying Lemma 14.51 again with d = d^ and ip = d* , we obtain 

\Eidoo) - E{d*)\ l - e < || - Ad* + f(d*)\\ = 0, (5.27) 
which together with (|5.25p yields ()5.12p . The proof is complete. □ □ 
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Remark 5.5. We note that in the assumptions 1 1 f o 1 1 h 1 < 1 an d \\do — cZ* ||h 2 — 1> ^ e bound 
1 is not essential and it can be replaced by any fixed positive constant M. In this case those 
constants in the proof of Theorem \5.S\ may also depend on M . □ 

Corollary 5.1 (Nonlinear stability). Suppose that n = 3 and the conditions in Case I are 
satisfied. Let d* E Hp be a local minimizer of E(d). Then d* is Lyapunov stable. 

Proof. We observe that in the proof of Theorem 15.21 oj can be an arbitrarily small positive 
constant satisfying oj < ^ min{<7, /?}, by our choice of o~\,o~%, we actually have shown that the 
local minimizer d* is Lyapunov stable. □ □ 

Remark 5.6. We can see from (|5.26p and (|5,27p that the asymptotic limit doo obtained in 
Theorem \5.2\ (ii) is also a local minimizer of E{d) (having the same energy level as d* ). Moreover, 
if d* is an isolated local minimizer, then d^ = d* and d* is asymptotically stable. □ 



6 Parodi's Relation and Linear Stability 

In Section 5.2 we have shown that Parodi's relation can be viewed as a sufficient condition for 
the nonlinear stability of the Ericksen-Leslie system ([l.lip — (|1.13p . It is still an open problem 
whether similar result holds true for the original Ericksen-Leslie system (jl.ip - (|1.3p . Alterna- 
tively, in this section we shall make a preliminary study to discuss the connection between 
Parodi's relation and linear stability of the original Ericksen-Leslie system (|l.ip - (|1.3p . 
For the sake of simplicity, we assume that 

p=l, fn=0, F = G = 

and the Oseen-Frank energy density function takes the simple form 

W = \\Vd\ 2 . 

Besides, we are interested in the bulk properties of the Ericksen-Leslie system and consider 
the problem in the whole space HL 3 , neglecting the boundary effects. We thus simply set the 
Lagrangian multiplier j3 = since /3 does not enter into the local field equations and can be 
determined through the boundary conditions, if any, on the director stress (cf. [8]). The nematics 
usually adopt a constant orientation in uniform shear flow. The analysis in [26} Section 6] shows 
that for a material that aligns in shear flow the viscous coefficients must satisfy 

|A*5 - Me| > |M2 - Mai- 
Here, we assume that (|1.16p and ()2.3p are satisfied, then we have 

|M5 - Mel > M3 - M2 > 0. (6.1) 

Consider the basic uniformly-oriented equilibrium state in which the material is at rest (zero 
velocity), the orientation is uniformly parallel to a constant unit vector n = (m, ri2, n^) T , the 
hydrostatic pressure p is a constant and the director tension (Lagrangian multiplier) 7 is zero. 
The equilibrium state is disturbed by perturbations with a small amplitude: velocity field v, 
director d + n, pressure p + p and director tension j. After a direct computation, the linearized 
equations of the Ericksen-Leslie system (|l.ip - (|1.3p for (y, d) are (cf. e.g., [7]) 

M3 + M6 _ M4 

0, (6.2) 



dvi H2 + M5 

-q£ + P,i - liiriinjn h niVj^i njn k v k:ij 

Ms - M2 „ ddij ddjj 



-njn k v i>k j - M2%-^ M3«-i 



dt ^ 1 dt 
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-A 



ddi 



diUi = 0. 



77ii - duj + 



Ai — Ao 



-rijVij 



Ai + A 2 



rijVjj = 0, 



(6.3) 
(6.4) 
(6.5) 



We study the behavior of infinitesimal, sinusoidal disturbances by a linear stability analysis. 
For this purpose, we seek plane wave solutions to the linearized system ()6.2|) — f)6.5[) of the following 
form (see e.g., [9]) 



d 

v 

7 
P 



ae 



l^gV - ^-{rav-x— lot) 
(JgV~ l(mu-x— u)t) 
■0^\/—l(m,u-x—u)t) 



(6.6) 
(6.7) 
(6.8) 
(6.9) 



where m is the complex wave number, v = (^1,^2,^3)^ is a given unit vector specifying the 
direction of propagation of the wave and oj is the complex frequency number, a, b are two 
constant vectors and C and D are two constants. Due to the constraint on the unit length of 
the director (|6.5p and the incompressibility condition (|6,3p . the constant vectors a and b satisfy 



n ■ a = 0, v ■ h = 0. 



(6.10) 



Let < 9 < ^ be a constant angle such that sin 9 = v ■ n. We consider the in-plane mode and 
deduce from (lUTUjl that (cf. [§]) 



a = A{y — n sin 9), 
b = B(n — z/sin#), 



(6.11) 
(6.12) 



where A and B are two constants. Inserting (|6.11|) and (|6.12|) into the linearized system (|6.2p - 
(|6.5p . after direct but tedious computations, we obtain 



(m 2 + y^Aiw) A 



-lm 



mujp(9)A 



m 2 g{9) 



q{0)B = 0, 
Au j B = 0, 



C v— lm + sin ( 



-T ^2 + 2 2 

-lwH m cos 1 



»' 2 (/'•-> - M2) D .... s 



L> + sin 6> m z + V-lAiw A + 



i? sin — fi2 m wA sin = 0, 
A2 — Ai 



lml? sin 9 = 0, 



where 



9(0) 



and 



2^i cos 2 sin 2 9 + (/i3 + /^6) cos 2 9 + m + (^5 — ^2) sin 2 
p(#) = H2 sin 2 9 — [13 cos 2 0, 

(Ai + A 2 ) cos 2 9 + (Ai - A 2 ) sin 

(Ai2 - ^3 + - Z^) COS 2 6> + (/X 2 - /i3 - ^5 + A*6) SU1 2 6 



(6.13) 
(6.14) 

(6.15) 
(6.16) 

(6.17) 
(6.18) 

(6.19) 
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Lemma 6.1. Suppose that (|1.16p . (16. ip are satisfied and M2M3 > 0. There exists a unique real 
solution 6>o S [0, f ] to the equations 

'pW = ' (6 20) 

j/ and onfo/ if Parodi 's relation (|1.17p holds. 

Proof. If (|1.17p holds, we have g(#) = 2p(9). Then it follows from (|6.ip and M2M3 > that the 
unique solution to (|6.20p is given by 

#0 = arctan ( */— ) G [0, — ]. 



Conversely, suppose 9q £ [0, f ] is the solution to (|6.ip . We discuss three subcases. 

Case 1: 9q = 0. It follows from (|6.18p that ^3 = 0. Then we infer from (|6. 19|) that 

M2 + M5 - M6 = and as a result, fj, 2 + M3 = M6 - Ms- 
Case #0 = §• I n this case we have M2 = and the proof is similar to Case 1. 
Case 3: < #0 < § . In this case it is easy to see that \ii 7^ 0, fi^ ^ 0. Since sin6*o 7^ 0, 

cos 6»o / 0, we deduce from ffTTKI) . (IBTTOj) that 

(Ai - A 2 )/u 3 = -(Ai + A 2 )/U 2 , 

which combined with (|1.16p yields (|1.17p . The proof is complete. □ □ 

In the remaining part of this section, we always suppose that ()1.16|) is valid. We make the 
following assumptions on the Leslie coefficients fi2, M3> M5> M6 : 

/i 6 > 0, fi 2 > 0, (6.21) 
fi 5 < min{/i 2 ,/i6}, (6-22) 
M3 = M6 - M5 + M2 - e, (6.23) 

where 

2(^6 - Ms)(M2 - Ms) 

u e <. mm < /Lt 6 - /Lt 5 , z/i 2 , 

Then we have 



< e < min <U 6 ~ Ms, 2/^2, V , _7 o h ( 6 - 24 ) 
[ 4/i 6 - 3/x 5 + 3^2 



Lemma 6.2. Under the assumptions ()6.21|) — ()6.24j) . f/ie Leslie coefficients satisfy conditions 
(|2.3|) and (|6.1|) . au£ Parodi' s relation (|1.17|) does no£ ZioZd. Moreover, there exists a unique 
solution #0 £ (0, § ) suc/i i/iai p(0o) 7^ and q(9o) = 0. 

Proo/. It easily follows from (fLT6D . (l6T23|) and (IBT241 that ([13]) is satisfied. Besides, $6?23\\ and 
(pT24"|) also imply that 

M6 ~ M5 < M3 ~ M2 + 2/i 2 = M2 + M3, 

so Parodi's relation (|1.17|) is not valid in this case. (|6.ip can be deduced from (|6.22p . (|6.24|) and 
(|2.3p in the sense that 

I Ms - M6| = M6 ~ Ms > M3 - M2 > 0. 

Finally, (j2.3|) and ()6.2ip yield that M2M3 > 0. Therefore, we can deduce from Lemma IBTTl that 
there exists an angle 



/ Me - Ms + M3 - M2 
arctan A / 

M6 - Ms - M3 + M2 



a ,,,an A / ^ e ( 0l E) (6 . 25) 
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such that 

p(e ) / and q(9 ) = 0. (6.26) 
The proof is complete. □ □ 

We further assume that 

< mi < i(2A*6-A*5 + M2), ( 6 - 27 ) 
0<^ 4 < ^(2^6-^5 + ^) cos 2 (6.28) 

where 6q is defined by ()6.25p . Then we can state the main result of this section: 

Theorem 6.1. Suppose that the Leslie coefficients (jbi, (mq satisfy the assumptions (|6.21|) - 
(|6.24|) , ()6.27p and (|6.28[) . Then the linearized Ericksen- Leslie system (|6.2p - ()6.5p admits unstable 
plane wave solutions. 

Proof. Let 9q be the angle obtained in Lemma 16.21 (cf. (|6.25p ). Taking 9 = 9q in the equations 
(|6.13p and (|6.14j) and using (|6.26p , we obtain that 



m 2 + V^IAicj) A = 0, (6.29) 

) 



mup(8 )A + ( m y o) - v^lu;^ B = 0. (6.30) 



Choosing 

A = 0, B = 1 and m 6 R, m / 0, 



we deduce from (|6.30p that 



u = -V^l (6.31) 

The (imaginary) constants C and Z? are determined by (|6.15|) and ()6.16p . respectively: 

C = m S ^ — [5(6*0) + ^2 - + ^5 cos 26*o] V-i, 

n A 2 - Ai r — . 

L> = v — lmsmfn. 

2 u 

(pHl]) implies that 

tan 2 9 = 2(/i6 = ^ " 6 > 2(2/*-^+/^ ^ 
Consequently, we deduce from (|6.27p . (|6.28p and (|6.32p that 

g(0o) 

= 2/xi cos 2 6*o sin 2 # + (^3 + n§) cos 2 9 + /M 

+(^5 - M2) sin 2 6> 
= 2^i cos 2 6*o sin 2 #0 + /M + (2^6 - ^5 + ^2 - e) cos 2 9 

+(^5 - /U2)sin 2 6> 



< cos 2 9 (2//i + /U 4 sec 2 6*o + [(2/u 6 - fi 5 + /i 2 ) + (a*5 - M2) tan 2 6*0] ) 

(2^ 6 - + /i 2 ) , 

< COS 9n 

2 u 

< 0. (6.33) 
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Thus, we obtain the following plane wave solutions (d,v,j,p) 

d = 0, 



(n — ^sin6*o)e x 



7 



-lmu-x- 



-lmu-x- 



which are unstable since g(0o) < 0. Here, we note that v — nsin#o 7^ due to the fact 6q G (0, ^). 
The proof is complete. □ □ 

Remark 6.1. Theorem \6.1\ indicates that for the nematic liquid crystal flow, if Parodi's relation 
(|1.17p does not hold, the original Ericksen-Leslie system (|l.ip - (|1.3p system can be (linearly) 
unstable. □ 



7 Appendices 

In this section we provide some detailed computations used in the previous sections. 



7.1 Least action principle 

The action functional takes the form 

fl 



A(x) 



1 



x t (X,t) 



\¥- T V x Ed (X)\ 2 + T(Ed (X)) 



JdXdt, 



where = Q is the original domain occupied by the material, E is the deformation tensor 
satisfying (|3.2p and the Jacobian J = detF = 1. The above expression includes all the kinematic 
transport property of the molecular director d. With different kinematic transport relations, we 
will obtain different action functionals, even though the energies may have the same expression 
in the Eulerian coordinate. 

We take any one-parameter family of volume preserving flow map 



x e (X,t) with x° = x, 



e=0 



and the volume-preserving constraint V x . • y = (or J € = detF e = 1). Applying the least action 
principle, we have 

dA(x £ 



such that 



T 



xt ■ ytdXdt 



dc 







e=0 



J Co 
T 



(¥~ T V x Ed ) 



'0 Jn 

' [ /(Mo) 
io Jn 

h + h+h, 



d 



(V x ed(x e ,t)) 



e=0 



dXdt 



dE e 
~de 



t=0 



d | dXdt 



(7.1) 
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where E e = E(x £ (X, t), t). Pushing forward to the Eulerian coordinate, we have 

r-T 



h 



xu ■ ydXdt 



o JQo 
T 



T 

o Jn t 



v ■ ydxdt 



o Jn 



(vt + v ■ Vv) ■ ydxdt, 



where is the domain occupied by the material at time t 
From the definition of E e , we have 

dE 



de 



do 



which implies that 
h 



e=0 



T 



I( Vy - V T y) - A( Vy + V T y) ) Ed,,. 



o Jn 

T 



(¥- T V x Ed ) : 



(F- T v x m ) 



d(F e 



o Ja 

T 



de 



¥- T V x 



V x Ed dXdt 



6=0 

dW 
de 



e=0 



dXdt 



[ f Vd : (-V T yVd) dxdt 
Jo Jn t 



T 



o Jn t 

T 



Vd: V 



Vy-V T y X 2 Vy + V T y 



o Jn 



2 Ai 
[V • (Vd Vd)) ■ ydxdt 



dxdt 



2 V Ai 



T 



[V • (Ad <g> d)] ■ ydxdt 
-\ (l + £ [V • (d ® Ad)] • ydxdt, 



Jflt 
T 



and 



Jtt 



l -(Vy-V T y)-^(Vy + V T y)) < ! 



dXdt 



T 

o Jn t 



A 2 



V • (/(d) d) 



+ 



o Jn 



i(l + ^)V-(d®/(d)) 



y dxdt 
• y dxdt. 



Inserting (j7.2|) . (j7.4j) and (|7.5p into (|7.ip . we arrive at 

/ / [v t + v-Vv + V ■ (Vd Vd) - V ■ a] ■ ydxdt = 
Jo Jn t 

A 2 



where 



1 



1 



(Ad - /(d)) d + i (l + ^)d (Ad - /(d)). 

2 V Ai / 



(7.2) 



(7.3) 



(7.4) 



(7.5) 



(7.6) 



(7.7) 



2V~ Ai- 

Since y is an arbitrary divergence free vector field, we formally derive the momentum equation 
(Hamiltonian/conservative part) after integration by parts 

v t + v ■ Vv = —VP - V • (Vd Vd) + V • <7, (7.8) 

where the pressure P serves as a Lagrangian multiplier for the incompressibility of the fluid. 
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7.2 Maximum dissipation principle 



Using the maximum dissipation principle [37-39], we perform a variation on the dissipation 
functional (half of the total rate of energy dissipation D (|3.7p ) with respect to the velocity v 
in Eulerian coordinates. If 5 V [^T>) is set to zero, we will get a weak variational form of the 
dissipative force balance law equivalent to conservation of momentum. Let v e = v + eu, where 
u is an arbitrary regular function with V • u = 0. Then we have 

1 dV(v e )_ 

6 = 







1 



-V 



de 



2 



\7v : Vudx + ii\ J d^Akpdpdi 



^^djdx 



d t + v -Vd - ftd+ —A d 
Al 



u ■ Vd 



Vu - V T u , A 2 Vu + V T u 



-d + 



A i 



d )dx 



+ [fJ-5 + Me + 



(A 2 ) 2 
Ai 



Q 



Aijdj 



d k dx 



:= h+h + h + h- 
Using integration by parts, we get 



M4 



(Au,u), 



h 



Using the transport equation (|3.6p of d and the incompressibility of u, we infer that 



m ( V • [d T ,4d(d <g> d)] , it 
~(a«5 + M6 + ^~~~) y (u k Vi(d k A^dj) + UiV k (Aijdjd k ))dx 
V • (d <g> Ad) , m) + (V • (Ad (g> d) , u 



1 / (A 2 ) 2 



Ad-f(d),u-Vd-±(l-^ 



V T ud 



u, -VF(d) + V • (Vd Vd) - V 



|Vdf 



1 + 



A 2 _ 
i 

A 2 



+ (l--£)(„,V-[(Ad-/(d))®d] 



A 



«,V- [d0 (Ad -/(d))] 

(«,V • (Vd0 Vd)) - ^ 2 (n,V • (N <g> d)) - fi 3 (u,V ■ (d(g>N)) 
-7/5 (u, V • (Ad<g)d)) - 7/6 (it, V • (dig) Ad)), 



with the coefficients 



(7.9) 



M2 



'/5 



2(Ai"A 2 ), 



A, 



(A 2 ) 2 



A 



M3 
7/6 : 



-2(Ai + A 2 ) 



A 2 + 



(A 2 ) 2 



It follows from the above calculations that 

1 dV 

2 de e=o 







V- (VdQVd)) - (u,V-<t). 



(7.10) 



(7.11) 
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The stress tensor a is given by 

a = [ii(d T Ad)d ® d + [i 2 N <3 d + [i 3 d ® N + + ft 5 Ad <3 d + fi e d <3 Ad, 
with constants 

^2 = ^(Ai-A 2 ), ^3 = -^(Ai + A 2 ), 

h = ^(A2 + + Me), Me = ^( _A 2 + /^5 + Me)- 
Since u is an arbitrary function with V-u = 0, we arrive at the dissipative force balance equation 



where the pressure P serves as a Lagrangian multiplier for the incompressibility of the fluid. 
7.3 Computation on the time derivative of A(t) 

Using (|l.lip - ()1.13p and integration by parts, due to the periodic boundary conditions, we have 



= -VP - V • (Vd Vd) + V • a, 



(7.12) 




-(Av,v t ) + (Ad-f,Ad t -f'(d)d t ) 



(Av,v Vu) + (Av, VdAd) - (V • a, Av) + — ||V(Ad - /)|| 2 



- (Ad - f,A(v Vd)) + (Ad - /, A(J2d)) - -p (Ad - /, A(Ad)) 




(7.13) 



First, we expand the third term on the right-hand side of ()7.13p : 



(V • a, Av) 




Using integration by parts and the fact that f2 is antisymmetric, we have 
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-fii / AkpdkdpViididjjViAijdx. (7.14) 
JQ 

By the incompressibility condition V • v = 0, we see that 

■/ o Jo 



-\\Avf. (7.15) 



Meanwhile, 



-M2 / Vi(djNi)ViVjVidx - fi 3 / V i(diN j)V iV jVidx 
JQ JQ 

fi2 / djNiA(Aij + Qij)dx + ^ 3 / diNjA(Aij + Qij)dx 
JQ JQ 

(^2 + ^3) / djiViAAydz + (M2 - V3)(N,Attd), (7.16) 



and 



-/i 5 / Vi(djd k A ki )ViVjVidx - n 6 / V i(did k A kj )V t V jVidx 
JQ JQ 

= djd k A ki A(Aij + ttij)dx + /i 6 / djd k A ki A(Aij - Q^dx 

JQ JQ 

= (Ms + / djd k A ki AAijdx + (/z 5 - /i 6 ) / djd k A ki AVL i: jdx 
JQ JQ 

= -(M5 + Me) / djd k V iA ki V lAijdx - (fj, 5 + ^ 6 ) / Vidjd k A ki ViAijdx 
JQ JQ 

-(^5 + Me) / djVid k A k iViAijdx + (fi 5 - ^(Ad, AQd) 
JQ 

= -(Ms + Me) / |(ijV^ji| 2 dx - (/i 5 + // 6 ) / V idjd k A ki V [Aijdx 
JQ JQ 

-(^5 + /^)/ djV L d k A ki V lAijdx + (n 5 - n 6 )(Ad,Arid). (7.17) 
•/ Q 

Next, using the d equation (|1.13p . we have 

(Ad-/,A(Q d)) 

= (Ad - /, AQ d) + 2 / (Ad* - f i )V l Q, ij Vid j dx 

JQ 

+(Ad-f,UAd) 
= -Ai / djNiAQijdx- X 2 (Ad, AQd) 

JQ 

+2 / (Adj - fi)Vin tj Vidjdx + (Ad - /, QAd) 
JQ 

= -Ai(iV, AOd) - \ 2 {Ad,And) - [ Vi{Adi - f^Vidjdx 

JQ 

+ / (Ad; - fijViQijVidjdx, (7.18) 

JQ 
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~(Ad — f,A(Ad)) 
\ AC AdS\ + i 

Ai 



A 2 (JV, A{Ad)) + (Ad, A{Ad)) 



A 



/ NiAAijdjdx + 2A 2 / NtfiAijVidjdx + \ 2 {N, AAd) 
JQ JQ 

I \Vl(Aijdj)\ 2 dx. (7.19) 
X l JQ 

Special Cancellations, (i) Due to (|1.16p . the first term on the right-hand side of (|7.18p cancels 
with the second term of the right-hand side of (|7.16p and the second term on the right-hand 
side of (|7.18p cancels with the fourth term of the right-hand side of (|7.17p . (ii) By Parodi's 
relation (|1.17p . the first term of the right-hand side of (|7.19p cancels with the first term of the 
right-hand side of (|7.16p . 

Concerning the fifth term on the right-hand side of (|7.13p . using the incompressibility of v, 
the fact Vd ■ f(d) = 'VJ-(d) and integration by parts, we obtain 

-(Ad-f,A((vV)d)) 

= -(Ad - f, Av • Vd) - 2 I (Ad, - fi)ViVjViVjdidx -(Ad-f,v VAd) 

JQ 

= -(At>,VdAd) + 2 I Vj(Adi- fi)ViVjVididx-(Ad- f,v-Vf). 

JQ 



Hence, 



- (Ad -f,A(v Vd)) + (At?, VdAd) 

+ (Ad - /, f'{d) \^-{Ad - f) + v • Vd - n d + ^Ad 
V La i Ai J 

= ^J^f'(d)\Ad-f\ 2 dx- \Ad-f,f'(d)(nd-^Ad^ 

+2 / Vj(Adi - fi)ViVjVididx -(Ad-f,v Vf). (7.20) 
JQ 

Collecting the above calculations together, we conclude that (|4.11|) holds. 
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